Competition between superconductivity and the pseudogap phase in the t-J model 
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The t-J model in the large N limit (N denotes the number of spin components) yields a pseudogap 
phase in the underdoped region which is related to a d-wave charge density wave (d-CDW). We 
present results for the doping dependence of the superconducting and d-CDW order parameters as 
well as for collective excitations in the presence of these two order parameters. We argue that the 
electronic Raman spectrum with B\ g symmetry probes the amplitude fluctuations of the d-CDW 
at zero momentum. 

PACS numbers:74.72.-h, 71.10.Hf, 71.27.+a 
The superconducting transition temperature T c in the cuprates shows a maximum at optimal doping 5 — 5 C and 
a monotonic decrease towards lower or higher dopings. On the other hand, the excitation gap in the one-particle 
spectrum, often called pseudogap, increases monotonically with decreasing doping, remains finite above T c in the 
underdoped region, and has d-wave symmetry. Generalizing the t-J model from 2 to N spin components it has been 
shown that such a pseudogap arises in the large N limit of the t-J model due to an instability of the normal or the 
superconducting state with respect to a d-CDW (sometimes also called flux or bond-order wave). In the following we 
present results for the symmetry broken states, namely the doping dependence of the superconducting and the d-CDW 
order parameters and the dynamics of d-wave charge density fluctuations in the presence of these order parameters. 
The relevance of a d-CDW for high-T c superconductors has recently also been discussed in Ref. [||. 
In the case of the t-J model the d-CDW order parammeter $ is defined as 

$ = -2i J/N c £ 7(k) (4 CT 5k+Q CT >. (1) 

ko- 

J is the Heisenberg coupling, , c are creation and annihilation operators for electrons under the constraint that 
double occupancies of lattice sites are excluded, N c is the number of primitive cells, (...) denotes an expectation value, 
and Q is the wave vector of the d-CDW. y(k) is equal to (cos(k x ) — cos(k y ))/2. Keeping only the instantaneous 
contribution in the effective interaction || the order parameter A for d-wave superconductivity is 

A = 2(J-K)/Ac$>(k)(ck T 5-ki). (2) 

k 

V c is a repulsive nearest-neighbor Coulomb potential which is needed to stabilize the d-CDW with respect to phase 
separation Q. In the presence of the above two order parameters the operators (c^ *, c_k,j., c^+q p c_k-Q,i) are 
coupled leading to the following Green's function matrix pi 
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-A(k) uj + e(k) i$(k) 
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(3) 



e(k) is the one-particle energy, e(k) = —(St + aJ)(cos(k x ) + cos(k y )) — 2t' 5cos(k x )cos(k y ) — fi, with a = 
l/N c cos (9a;)/( e ( c l))- / is the Fermi function, S the doping away fom half-filling, fi a renormalized chemical poten- 
tial, t and if are nearest and second-nearest neighbor hopping amplitudes, to a (complex) frequency, and k = k — Q. 
<I>(k) and A(k) are equal to <E>7(k) and A7(k), respectively. 

Expressing the expectation values in the order parameters by G and using Eq.(^) one obtains coupled equations 
for the two order parameters. The thick solid and dashed lines in Fig. 1 show the numerically determined doping 
dependence of $ and A at zero temperature, calculated for t'/t = —0.35, J/t = 0.3, V c /t = 0.06, and Q = (it, 7t). 
The energy unit is t. The thin solid and dashed lines show $ and A for non-interacting order parameters putting the 
second order parameter to zero. In this case the onset of $ is much higher than in the interacting case, especially 
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FIG. 1. Order parameters <E> and A as a function of doping in units of t at T = in the interacting (thick lines) and 
non-interacting (thin lines) case. 



if t' is nonzero. This reduction in the onset of <E> reflects the competition of the two order parameters. The supercon- 
ducting order parameter A increases monotonically with decreasing doping in the non-interacting case but is heavily 
suppressed by $ below the onset of $ in the interacting case. For 0.05 < 5 < 0.15 both order parameters are non-zero 
so that superconductivity and d-CDW coexist throughout this region. 

The eigenvalues of the matrix G~ l give the one-particle energies in the presence of the two order parameters. The 
corresponding density of states p(uj) is shown in Fig. 2 for three different dopings. In the upper diagram, corresponding 
to the overdoped case with $ = 0, the usual density for a d-wave superconductor is seen. The lowest diagram in 
Fig. 2 corresponds to the underdoped region dominated by the d-CDW. The zero in p(u>) occurs somewhat above 
the chemical potential lo — and p(co) is rather asymmetric around this point. The diagram in the middle of Fig. 2 
describes a slightly underdoped case where both order parameters are of similar magnitude. p(oj) reflects here both 
gaps and possesses additional structures due to the geometry of the two-dimensional Fermi surface and the Umklapp 
processes induced by Q. 
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FIG. 2. One-particle density of states for three different dopings. 




The above calculations were performed assuming a commensurate Q = (tt,it). In reality the d-CDW is at low 
temperatures and a finite doping incommensurate [|IJ with four inequivalent wave vectors Q = ±(7r, q), ±(q, tt) with 
q different from n. Writing q = n — x the incommensuration x was approximately calculated from the leading 
instability in k-space using a linear approximation, x is shown in Fig. 3 as a function of S. The left diagram in this 
figure corresponds to A = 0, i.e., the thin line in Fig. 1. In this case x assumes the rather large value of about 
0.6 at the onset of the d-CDW, but then decays rapidly to about 0.1 with decreasing doping. The right diagram in 
Fig. 3 illustrates that superconductivity not only depresses the onset of the d-CDW but also the incommensuration x 
resulting in a practically commensurate d-CDW in the underdoped regime. In our calculations x has been determined 
from that momentum where the d-wave charge susceptibility in the superconducting or normal state shows the 
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strongest divergence. 
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FIG. 3. Deviation x of one wave vector component of the d-CDW from n as a function of doping in the absence (left diagram) 
and presence (right diagram) of superconductivity. 

Let us denote by x(q, w) the response function associated with the d-wave density operator p(q) = 
V^Eka^t^^+q^iff- At large N x is given by x(q, w ) = X (0) (<I> w )/(! + J X (0) (q, <*>)) which holds exactly at 
q = (0, 0) and (71", n) and in a good approximation for a general q. For a general q the explicit expression for is 
rather involved, so we give it here only for the special case q = (0, 0), 

X (0) H - fll.llH - ^2,2l(«) - Pw.SlH + Pl4,41 M, (4) 



PijM{u) = — ^7 2 (k)Gij(iw n + w, k)Gfcj(iw„,k). (5) 

C k,n 

The left diagram in Fig. 4 shows the negative imaginary part of x*- -* (dashed curves) and of % (solid curves) at 
q = (0, 0) and three different dopings. The curves illustrate that the Heisenberg term as the residual interaction 
strongly scatters the quasi-particle excitations across the pseudogap shifting most of the spectral weight down into 
bound states inside the gap. In the overdoped case (upper diagram in Fig. 4) this bound state describes an exciton 
state inside the superconducting gap. In the strongly underdoped case (lowest diagram) the bound state corresponds 
to the amplitude mode of the d-CDW probed at a 
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FIG. 4. Left diagram: Correlation function for d-wave density fluctuations at q = (0, 0) in the free (dashed lines, called B{1 ) 
and interacting (solid lines, called Bi a ) case; Right diagram: dispersion of the main peak in the d-wave density fluctuation 
spectrum along (q, q) for three different dopings. 
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wave vector (tt, tt) away from the wave vector of the d-CDW, i.e., at q = (0, 0). With decreasing doping the frequency 
of the peak increases monotonically and thus is tight to the total pseudogap. The solid curves in Fig. 3 are proportional 
to the Big spectra of electronic Raman scattering in the cuprates and are in good agreement with the experimental 
data In particular, one can understand why the B\ g peak does not probe the superconducting but the pseudogap 

i- 

In the right-hand diagram of Fig. 4 we plotted the position of the main peak in the d-wave fluctuation spectrum 
along q = (q,q) for three different dopings. In the overdoped (6 = 0.25) and the slightly overdoped (S = 0.15) cases 
only one well-pronounced peak in the spectrum was obtained describing collective d-wave density fluctuations. Its 
dispersion is rather weak away from (tt, tt). Near (tt, tt) a well-pronounced soft mode develops if the doping approaches 
the onset of the d-CDW from above. Below the onset the soft mode hardens with decreasing doping and its peak 
position is near or somewhat below the value 2$ in analogy with BCS-theory. In the strongly underdoped case with 
8 = 0.077 the momentum region around (tt, tt), where the peak position is roughtly given by 2$, is rather small. In 
the neighborhood of q ~ 2 the fluctuation spectrum develops two peaks and the spectral weight shifts from the upper 
to the lower peak with decreasing momentum. The lower peak shows only little dispersion towards smaller momenta 
and lies somewhat above the corresponding peaks at larger dopings. 

In conclusion, we found that the underdoped regime of the t-J model at large N is characterized by the competition 
between d-wave superconductivity and a d-wave CDW. This competition is especially strong for a finite t'\ For 
instance, for t' /t = —0.35 the d-CDW onset shifts from about 6 = 0.29 to S = 0.14 due to superconductivity. The 
wave vector of the d-CDW is very close to the commensurate value (tt,tt), except near its onset. The d-wave density 
fluctuation spectrum is determined by collective effects. It exhibits a pronounced soft mode behavior near the onset 
of the CDW but shows only little dispersion far away from (tt, tt). The fluctuation spectrum at zero momentum agrees 
well with data from electronic Raman scattering, in particular, with respect to its dependence on doping. 

The authors thank Secyt and the BMBF (Project ARG 99/007) for financial support and P. Horsch for a critical 
reading of the manuscript. 



[1] E. Cappelluti and R. Zeyher, Phys. Rev. B 59, 6475 (1999) 

[2] S. Chakravarty, R.B. Laughlin, D.K. Morr, and Ch. Nayak, Phys. Rev. B 63, 94503 (2001) 

[3] R. Zeyher and A. Greco, Eur. Phys. J B 6, 473 (1998) 

[4] C. Kendziora and A. Rosenberg, Phys. Rev. B 52, R9867 (1995) 

[5] R. Zeyher and A. Greco, to be published. 



4 



